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OPTIMAL TRANSPORT WITH COULOMB COST. 
APPROXIMATION AND DUALITY. 

LUIGI DE PASCALE 


Abstract. We revisit the duality theorem for multimarginal optimal trans¬ 
portation problems. In particular, we focus on the Coulomb cost. We use a 
discrete approximation to prove equality of the extremal values and some care¬ 
ful estimates of the approximating sequence to prove existence of maximizers for 
the dual problem (Kantorovich’s potentials). Finally we observe that the same 
strategy can be applied to a more general class of costs and that a classical 
results on the topic cannot be applied here. 


1. Introduction 


be a 


This paper deals with the following variational problem. Let p G 
probability measure (p will be called electronic density) and let 

c{xi,...,xn) = 

l<i<j<N 

be the Coulomb cost. Consider the set n(p) = {P G P(M^'^) | vrJP = p}, where 
TT* denotes the projection on the z—th copy of and ttJP is the push-forward 
measure. We aim to minimize 


Xi 


Xi 


min / c{xi,... ,XN)dP{xi,... ,Xn). (1.1) 

n(p) jRivd 

This problem fits in the general framework of multimarginal optimal transporta¬ 
tion problems. In particular it is the multimarginal optimal transportation prob¬ 
lem with all the N marginals coinciding with p and with Coulomb cost. 

In contrast with the classical two-marginal optimal transportation problems, the 
theory of multimarginal optimal transportation problems is still at the beginning 
and many relevant open problems need to be studied. Some general results are 
available in n HZl ESI ED E3], results for special costs are available, for example 
in [ID for the quadratic cost with some generalization in [ID, and in [5] for the 
determinant. Some new applications are emerging for example in |12j . 
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In the particular case of the Coulomb cost there are also many other questions 
related to the applications. Recent results on the topic are contained in |2l El El 
and some of them will be described better in subsequent sections. 

The literature quoted so far is not at all exaustive and we refer the reader to 
the bibliographies of the cited papers for a more detailed picture. However, in 
the author’s opinion, at the moment the wealth of problems obscures the body of 
known results. In this paper we will focus on the Kantorovich duality for problem 
fll.ip . We denote by the product measure p ® ■ ■ ■ ® p. We will prove that 

N—times 


min 

n(p) 



c{xi, . . . , XN)dP{xi, . . . , Xn) 


sup{W 



u E Lj, 


m(xi) H- \-u{xn) < c{xi,... ,xn), p^ 


a.e.} (1.2) 


and that the right-hand side of fll.2p admits a maximizer which is, in particular, 
bounded. Thanks to the symmetries of the problem we also have that the right- 
hand side of fll.2p coincide with 


N 

sup{^ / Ujdp : Ui E Lp, Ui{xi) + ■ ■ ■ + un{xn) < c{xi,... ,xn), p^ — a.e.}. 

1=1 d 


Infact, this last sup is a priori bigger then fll.2p . Since for any admissible N- 


N 


tuple {ui, ... ,uisf) the function u{x) = — Ui{x) is admissible for the previous 


i=l 


problem, equality holds. 

Some of the basic ideas originated in the paper El- However, in Remark 13.141 
we will show that Theorem 2.21 of na does not apply to the Coulomb cost in 
order to prove that a maximizer of the dual problem exists. The tools will be 
T-convergence and some careful estimates of the maximizers of the approximating 
problems. 

We remark that a necessary but not sufficient assumption for fll.ip to be finite 
is that p is not concentrated on a set of cardinality < N — 1. We will assume this 
whenever needed. 

We will adopt the notations x = (xi,..., xtv) G so that Xj E for each 
z G {1,..., N}. And p—inf will denote both the essential p infimum or the essential 
p^~^ infimum depending on the number of variables involved. 

1.1. Motivations. The main object in the quantum mechanical modeling of a 
particle with N electrons is a wave-function tjj, i.e. an element of 


A:={il;EH\{R^xZ2f,C) : II^IU^ = !}■ 


The space (M^ x Z 2 )^ is the configuration space for the N electrons. In fact, 
the state of each electron is individuated by the position in and the spin in 
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Z 2 . While it is easy to understand the interpretation of the position variable, the 
notion of spin is slightly more difficult to grasp. We may roughly explain it as 
follows; when a magnetic held is applied to the electron of position tq velocity uq, 
the electron may be dehected in two different ways which are characterized by the 
two possible spins. The name is then reminiscent of the behavior of a spinning 
top which hits a wall with a certain translational velocity in a point po and may 
rebound in two different directions depending on the whirling direction. However 
this analogy, although suggestive, is not the historical one nor can be considered 
physically meaningful. 

It is easy, once we free ourself from this last analogy, to imagine the possibility 
of a spin variable with values in or other spaces. 

The quantity 

|i/>((ti,q;i), ..., (a;jv,a7v))P 

represents the probability that the N electrons occupy the state ((ti, ai),..., 
(ttvjCItv)) and then, since the electrons are indistinguishable tf) satishes 

1(^cr(i) 9 ) 9 * * * 9 (^f7(7V) 9 ))| |i/^((ti,o;i),..., , Dtv) ) I 

for all permutations a of the N electrons. This brings to introduce the distinction 
between Fermionic and Bosonic particles, however such distinction will not have 
relevance here since we will discuss duality for a relaxed problem. The distinction 
may be instead relevant when discussing the relaxation process. A particle is 

Fermionic if 


l/^((^cr(l)9 ^Cr(l))9 * * * 9 (^Cr(A^) 9 ^f7(7V))) (^1 9 ^ 1)9 • • * 9 (^Y9 ^ v) ) 9 

and Bosonic if 


^((^f7(l)9 ^cr(l))9 * * * 9 (^f7(7V) 9 ^Cr(A^))) ^((^l9 ^ 1)9 * * * 9 (^Y9 ^y))* 

Since the spin will not play any role in the rest of the paper we will consider 
wave-functions depending only on the positions of the electrons. Including the 
dependence on the spin is just a matter of adding a summation over the two 
possible values in the formulas (below we give an example in the case of kinetic 
energy). In the simplest situation the electrons move with a certain velocity while 
interacting with the nuclei (or, equivalently, with an external potential) and with 
each other. The interaction with the nuclei and the interaction between electrons 
are of Coulombian nature. So if we assume that there are M nuclei with charge 
Zi and position Ri the interaction potential in the position x will be v{x) = 

Then the energy of a particular state of the particle is made of 

three parts: 


-E 


Z, 


X 
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the Kinetic ener 




/M3JV 


the electron-nuclei interaction energy 


N 


KeM = ^ / v{Xi)\'llj{xi, . . .XN)\‘^dXi . . .dXN, 


2=1 


/R3iV 


the electron-electron interaction energy 

1 


Veebl^] = [ 


l<i<j<A ' 


Xj — Xj 


, XAr)phXi . . - dXM- 


A relevant quantity in quantum mechanics, quantum chemistry and other dis¬ 
ciplines is the ground state energy of an atom 

Eq = m.mT[ip] -f Veeii’] + KeM- 


Remark 1.1. The special form of the electron-nuclei interaction is widely irrelevant 
in the expression of Vee above. So one could just include a term of the same form 
with an external potential v. 


In order to compute (numerically) Eq one needs to solve the Schrodinger equa¬ 
tion in ?)N dimensions and this is very costly even for a small number of electrons. 
A less costly alternative is represented by the Density-Functional Theory in¬ 
troduced hrst by Hohenberg and Kohn in [T6j and then Kohn and Sham in |18| . 
At the beginning of the theory the mathematical foundations of DFT were very 
weak. The paper which started to put DFT on solid mathematical foundations is, 
to my knowledge, [T9] . 

The main object in DFT is the electronic density p which is obtained by inte¬ 
grating out — 1 particles 

p{x) = / |'i/'(x, X 2 , . . . , X7v)pdX2 . . . dXTv 

JK3(JV-1) 

and it represents the probability distribution of a single electron. The relation 
between ip and p will be denoted hy ip p p. In particular p is always a probability 
density on independently of the number of electrons. The basic idea of DFT 
is to express everything in term of the electronic density p instead of the wave 


1 


if we want to include the "spin coordinate" in the state of the particle then 


T[iP] = 



|V'b((a:i,ai), ..., {xN,aN))\^dxi . ..dxN- 
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function ?/;. It is easy to express the electron-nuclei interaction energy in terms 
of p. We have 

KeM = N [ v{x)dp{x). 

JR3 

Then, following Hohenberg and Kohn, we write 

Eq = mm{FHK{p) + N [ v{x)dp{x)} (1.3) 

P 

where 

Fhk{p) = min{T[i/;] + WeM}, 

Pip 

is called the universal Hohenberg-Kohn functional [universal since it does not de¬ 
pends on the specihc particle). At the beginning of the DFT several mathematical 
questions needed to be settled in order to have a meaningful mathematical theory. 
For example: Which is the correct minimization domain? Does the functional 
Fhk enjoy some properties which make the variational problem amenable? And 
so on. As we said, Lieb in [19] started to look at these questions. In particular, 
Lieb described the set of admissible p 

H = {p\0<p, j p=l, 

Writing explicitly Fhk is not possible, then Kohn and Sham considered 

Fks[p) = mm{T[i;]}. 

Pip 

In the Bosonic case it is possible to express 

Fks{p) = N j \V^p\^dx. 

In the fermionic case 

j IVVpP* < < (4ir)=]V" J\VXp\^dx 

(But the optimal constant in the inequality above, to my knowledge is still not 
known). Then one may rewrite 

Fhk = Fks + F^c 

where F^c is called exchange-correlation energy it needs to be determined or at 
least approximated and it is the term which keeps into account the interaction 
between electrons. 

Let us denote by €.[p) the minimal value in fll.ip . Since for every if) \r p the 
measure \'ip\^dxi... dxn G n(p) we have 

€{p) < KM- 
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The functional £(p) appear in DFT in several ways. The hrst and perhaps the 
most elementary appearance is in the estimate from below 

Fks{p) + ^{p) < Fhk{p)- 

This estimate is useful in the variational study of DFT and is is also the basis 
for the so called KS-SCE DFT (Kohn and Sham, Strictly Correlated Electrons 
Density-Functional Theory) in which €.{p) is interpreted as an approximation of 
Fxc for particles in which the electron — electron interaction is relevant. This is 
a very rapidly developing domain started in |24| (see also [H [261 [25] but we are 
far from being exhaustive). 

Again C(p) appear when considering the SCE-DFT (Strictly Correlated Electrons- 
Density-Functional Theory), which is the analogous of DFT in a regime in which 
the electron — electron interaction is preponderant. This means writing 

Fhk = + Fkd, 

where Fkd is called kinetic-decorrelation energy, it needs to be determined and it 
is the term which keeps into account the kinetic energy of the particle (see, for 
example, [TTl [13] ). 

The third appearance is given in [S] where the authors proved that C is the 
semiclassical limit of the DFT in the case of a two electrons particle. It is not 
clear if this last result also holds for more than two electrons. 

Finally we quote a direct application of the duality theorem we are consid¬ 
ering. In the numerical approximations of all the problems above one need to 
compute €{p). Before the link with Optimal Transport problems was discovered 
this approximation was done by computing the co-motion functions which are the 
analogue of the optimal transport maps. However, more recently, Mendl and Lin 
in [20] introduced a numerical method which use the Kantorovich potential to 
compute the value. 


2. Preliminary results and tools 

2.1. Definition of F-convergence and basic results. A crucial tool that we 
will use throughout this paper is F-convergence. All the details can be found, for 
instance, in Braides’s book [2] or in the classical book by Dal Maso j^. In what 
follows, {X, d) is a metric space or a topological space equipped with a convergence. 

Definition 2.1. Let (Fn)n be a sequence of functions X i— )■ M. We say that (Fk)n 

r 

F-converges to F and we write Fk —)■ F if for any x E X we have 

n 

• for any sequence {Xn)n of -A converging to a: 

lim inf F„(a;„) > F{x) (F-liminf inequality); 
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• there exists a sequence (xn)n converging to x and such that 
limsup< F(x) (F-limsup inequality). 

n 

This dehnition is actually equivalent to the following equalities for any x G X: 

F(.) ^ i„f {.in, i„f ^ ^ ,nf {.nn.,p ^ h 

The function x i—)■ inf |liminf Fn{xn) : —)■ x| is called T-liminf of the sequence 

{Fn)n and the other one its T-limsup. A useful result is the following (which for 
instance implies that a constant sequence of functions does not T-converge to itself 
in general). 

Proposition 2.2. The T-liminf and the T-limsup of a sequence of functions {Fjfjn 
are both lower semi-continuous on X. 

The main interest of T-convergence resides in its consequences in terms of con¬ 
vergence of minima: 

Theorem 2.3. Let {Fn)n be a sequence of functions X —)■ M and assume that 

p 

Fn —)■ F. Assume moreover that there exists a compact and non-empty subset K 

n 

of X such that 

Mn G X, inf = inf 

X K 

is equi-mildly coercive on X). Then F admits a minimum on 
X and the sequence (inf Fn)n converges to minX. Moreover, if {xn)n 'Is a sequence 


having a limit x, then F{x) = inf F. 


(we say that {Fn 
X and the sequt 

of X such that 

\imFn{xn) = hm(inf F„ 

n - -ix 

and if (x 0 („))„ is a subsequence of (x^ 

Since we will deal also with dual problems we need the analogue of T-convergence 
for maxima, which is sometimes called r’''-convergence. 

Definition 2.4. Let (Fn)n be a sequence of functions X ha M. We say that {Fn)n 

I p + 

T^-converges to F, and we write F^ —> F if for any x G X we have 

n 

• for any sequence (xn)n of X converging to x 

limsupF„(x„) < F(x) (r’''-limsup inequality); 


• there exists a sequence (x„)„ converging to x and such that 
liminfF„(x„) > F(x) (T’^-liminf inequality). 

n 

The natural generalizations of Prop. 12.21 and Th. 12.31 hold with upper semi¬ 
continuity, maximum values and maximum points replacing lower semicontinuity, 
minimum values and minimum points. 
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2.2. Permutation invariant costs. The Coulomb cost c as well as the approx¬ 
imations that we will consider are permutation invariant in the sense that 

I^(^cr(l) 7 • • • 7 ^a(N) ) , . . . , ^7v) 7 G (S 

Several simplihcations are permitted by this invariance. 


3. Duality 


Denote by 


. X _ 1 1 a X E K, 
Xk\x) < Q otherwise. 


Introduce the set of elementary functions 

k 

= {v? : —)■ M measurable : if{x) = CiiXA^x---xA^} 

i=l 

for suitables constants Uj and Borel sets A], C M^. We consider the following 
approximation of c. 

Proposition 3.1. There exists a sequence of costs Cn such that 

(1) 0 < < c, 

(2) G 

(3) Cn Z ' C, 

(4) Cn is I.S.C., 

(5) Cn is permutation invariant. 

Proof. Let (ai,..., aZ G and dehne a — 1 := {k — l,p — 1, q — 1) for a = 
{k,p,q) G Z^ and (with a slight abuse of notations) {x,y) := {xi,yi) x (x 2 ,i/ 2 ) x 
( 2^37 Vs) foi' x,y E Then Vn G N dehne the interior of the diadyc cube by 

, N ,ai —1 ai, .ttAT —1 utvn 

Q iai ,..., otv) = (- 7 —) X • • • X (-, —). 

V 17 1 1'! j \ 2n ’ 2’’^' ' 2” * 2 ”"' 


and for all positives R 


Dehne c„ as follows 


Q,{R) = [-R, R] X ■ ■ ■ X [-R, R] 


3A-times 


inf c(x) if X G Q^{ai, ..., aZ and ..., aZ C Q(n), 

Cn(x) = -( Q"(aiv,«iv) 

0 if X ^ Q(n)- 

Then extend c„ to the entire space by relaxation i.e. 

Cn(x) = inf{hminf c(x^) : lim x^ = x, and 


k^oo 


k^oo 
.k 


yk, x^ E Q^{ai, ..., On) for some (ai,..., on) E Z^^}. 



DUALITY FOR THE COULOMB COST 


With this dehnition properties (1), (3) and (4) above are straightforward. To prove 
(2) it is enough to remark that is constant on i dimensional, relatively open 
faces of , oat) and that such faces are products of intervals (which may 

degenerate to a point). Finally (5) holds since 

inf c = inf c 

for every a G . □ 


Remark 3.2. Without loss of generality we can write 


Cn(^) 'y ^ ^iXAj (^l) ■ ■ ■ XaF (^y) 


2=1 
ifc ^ Ak 


with sets Af such that if n 7 ^ 0 then ylf = for all /c G {1 ,..., A^}. 


N 


Dehne iF(c„, •) ; ) —)■ [ 0 , +cxd] as 

K{cn, P) = I j ^ ^ 

I +CXD otherwise, 

and D{cn, •) : —)■ M as 

n) = I ^ j if H-1-n(xAr) < Cn(Ti, . . . ,XAr), p 

I —00 otherwise. 

Also in this case, maximizing D{cn,v) is equivalent to maximize 

N 

/ Vidp if ni(Ti) H-h vn{xn) < Cn{xi, ..., xn), - a.e. 

i=i 

—cx) otherwise. 

Proposition 3 . 3 . The functionals K{cn, ■) are equicoercive and 

K{c^,-)^K{c,-), 

with respect to the w*-convergence of measures. 


a.e. 


Proof. Equicoercivity follows from the fact that n(p) is zc*-compact. Since K{cn ,) 
is non-decreasing (in n) and K{c, ■) is l.s.c., F-convergence is a standard fact. We 
report the proof for the sake of completeness. Let Pn ^ P and hx m G N. For 
m < n 

K{Cm,Pn) < K{Cn,Pn)- 
Since Cm is lower semi-continuous 

K{cm,P) < liminf K{cm,Pn) < liminf K{cn,Pn)- 

n^oo n^oo 
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And since K{cm,, P) —^ A'(c, P) we obtain 

K{c,P) < \im ini K{cn,Pn)- 

n—^oo 

For what concerns the F — lim sup inequality, since c^, ^ c it is enough to consider 
Pn = P for all n G N and to apply the Monotone Convergence Theorem. □ 

Lemma 3.4. Assume that p is not concentrated on a set of cardinality < — 1. 

If u E Lp and 

u{xi) + • • • + u{xn) < c(xi,..., xn), — a.e., 
then there exists /c G M such that u < k p-a.e. 


Proof. We have 


u{x) < p — inf {c(x, X2,..., xat) — m(x 2) — • • • — n(xAr)}, p — a.e. 

X2y...,XN 


Since u E L^p and p has not < — 1 atoms, we may consider ... ,x^ with 

Hi 7^ Hj if ^ 7^ J such that u{xf) E M. We denote I = min{M(xi),..., m(xjv)}- 
Consider 0 < r such that 5(Xj, r) H B{Xj,r) = 0 if z ^ j. Then any x G can 
belong to B{x^,r) for at most one i. We may suppose without loss of generality 
that z = 1. It follows that 


m(x) < c(x,X2,... ,x^) — izfe) 


u{xn) < 




□ 


Remark 3.5. In the previous proof k clearly depends on u through the constant I 
and on p through the constant r. 

Proposition 3.6. The functionals 

D{cn, •) D{c, •), 

with respect to the weak-L^ convergence. 

Proof. We start with the F’*' — lim inf inequality. Let u E L^p such that —oo < 
D{c, u) then 

Zz(Xi) +-h Zz(xAr) < c(Xi, . . . , Xat), p^ — 0.6. 

and by Lemma 13.41 zz is bounded above. For an arbitrary 0 < £ there holds 

D(c, u) — e < N / udp. 

Jrs 

Moreover by Lusin’s theorem there exists X C compact such that u\k is con¬ 
tinuous, p(U.^ \ K) < e and 
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Since Cn c and are l.s.c we may apply a Dini’s type theorem and we obtain that 
there exists Uq such that if uq < n then 

u{xi)-\ - \-u{xn) - £ < Cn{xi,... ,xn), — a.c. (3.1) 

on K X ■ ■ ■ X K. We consider 

_ s 

u{x) := u{x) - ^ - N{supu)xm3\K{x). 

Since 0 < c„ and (13. Ij) holds we have, for no < n, 

u{xi) H-h u{xn) < Cn{xi ,..., Tat), - a.e. 

Moreover if M G N and N‘^ sup u < M then 


D{c,u) — {2+M)e < N / udp—e—N‘^esnpu<N / udp—e—N‘^p{M.^\K)snpu = 

Jr3 Jr3 

= N f udp. (3.2) 

Jr3 

For what concernes the F’*' — limsup inequality assume that ^ n in and 
assume without loss of generality that 

Un{xi) H-h Un{xN) < Cn{xi ,..., xn), P^ “ a.e. 


Then 

and 


Un{Xi) H-h Unix^) < c{Xi, ...,Xn) 


lim N 

n^oo 


Undp = N 



□ 


In the next subsection we will prove the needed compactness property. 

3.1. Estimates of the approximating Kantorovich potentials and conclu¬ 
sions. 

Lemma 3.7. For all n G N 

max Zl(cn, •) = min K(cn, •). (3.3) 

Moreover D{cn, ■) admits a maximizer (n",..., u^) 

Proof. The proof is revisited from im. 

kn 

Since we have written Cn(x) = OiXA^i^i)... Xa^{xn) we may hnd sets Xi,..., 

i=l 

each of + 1 elements, 

—)■ Xi 


and 


Cji '. X\ X • • • X Xjv —i K 
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such that 

Cn(x) = Cn{ipi{Xi), . . . , <Pn{xn))- 

If we dehne y? : —)■ Xi x ■ • ■ x Xj^j by 99 = and pi = ipi E V{Xi), 

then we have that for all 7 G n(/ 9 ) we may consider 7 := ip^'y G n(/ 9 i,..., p^) and 


Cnd-y = / Cn(v9(x))d7 = 


Cnd-y. 


IXiX---xXn 


We remark that since c„ and 7 ^ may be identihed with elements of ]^fcn + l X ... X 
I^fcn+i with an element of the problem may be reformulated as 

follows; 


mm c • 7 , 
Pil = Pi, 
0 < 7, 


(3.4) 


where the Pi form the projection matrix. Thus the problem is a linear program¬ 
ming problem whose minimum value coincides with the maximum value of the 
dual problem 


max p ■ u, 
P^u < c^. 


(3.5) 


It remains to identify u with a Wtuple {ui^... ,un) of elementary functions in 
T(M^^) which maximises 

N 

max{ E Uidp : UiEL^, Ui{xi)-\ - \-Un{xn) < Cn{xi,... ,Xn)}- (3.6) 

i=i d 


□ 


We now prove a first property of minimiser P of K{c, ■). This requires additional 
assumptions on p. 

Since the pointwise transportation cost diverges when two (3—dimensional) co¬ 
ordinates get too close it is useful to introduce the following notations: 

Da := {x = (xi,..., Xtv) : |xj — < a for some 

will be the closed strip around the diagonals, 

G ;= {x = (xi,..., xat) : x* 7 ^ Xj if f 7 ^ j}- 

For simplicity we pose Dq = D. Then G = \ D. Finally we introduce a 

notation for a "cube" with rounds 3—dimensional faces 

Q(x,r) := {y = (yi,... ,i/Ar) : max |yi - Xi| < r}, 

I 

we will omit the center x when it is the origin. 
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Proposition 3.8. Assume that p does not have atoms. Let P G n(/9) be a mini- 
mizer. Then, for all r > 0 there exists a{r) such that 

P{Da(r) n Q{r)) = 0 

We first prove an elementary lemma 


Lemma 3.9. Assume that p does not have atoms. Let P G Il{p)be a plan with 
finite cost and let x G sptP. Then there exist x^,..., x^ G sptP such that 

( 1 ) x^...,x^GG, 

(2) x] y^xl ifk-fii ora^j. 


Proof. First remark that if P has hnite cost then P{D) = 0 and then P-a.a. points 
belong to G. Then £x a vector a G and consider the set Xk = {x G : 

Xi = a}. By dehnition of marginals and since p does not have atoms 

P(Xi) = p({a}) = 0. 

Then, starting from x^ = x we may choose 


G sptP \ [D U U i=i,...,N X^k)). 


□ 


Proof, (of Proposition I3.8p Assume that x^ = (x},..., G DHsptP and choose 
points x^,...,x^ in sptP as in Lemma 13.91 For all choices of positive radii 

ri..., Tat 

P(g(x*,r,)) >0. 

And later on we will choose suitable rfs. Denote by P* = P|Q( 2 ,i and choose 
constants Aj G (0,1] such that 

= ■ ■ • = AatIPatI 

We then write 


P = AiPi + ■ ■ • + AatPat + Pr (Pr is the remainder). 
We estimate from below the cost of P as follows 

N 

C{P) = C{PR) + J2>^^C{P^)> 

i=l 

N N ^ 

C{Pr) + EuEE 

i=\ k=\ k<j 


XI 


x] I + 2rF 


I PI 


Consider now the marginals z/|of AjPj and build the new local plans 


Pi = I/} X ^2 X . . . P 2 = vl X 1^2 X . . . ulf, . . . , Prt = X X . . . V 


N-1 

N ■ 
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To write the estimates from below it is convenient to remark that we may also 
write: Pj = x ... x ... where we consider the lower index (mod N). 

Then consider the new transport plan 

P:=Pr + Pi + --- + Pn. 


It is straightforward to check that the marginals of P are the same as the marginals 
of P. Moreover \Pi\ = Xi\Pi\. So we can estimate the cost of P from above. 

N 


C(P) 


C(PH) + 5^C(Pj)< 

N N 

i=l k=l k<j ' 



The hnal step is to choose r* for z = 1,..., so that the 3—dimensional faces of 
the cubes Q(x*,ri) do not overlap and 


N N 

E(EE 

2=1 k=l k<j 


X 


k-\-i—l 


Xa 


\ -Tk- 


N N 

<E(EE 

2=1 k=l k<j 


XI 


x* I + 2ri 


This hnal condition contradicts the minimality of P and it is feasible because, 
since x] = x] for some i and j the right hand side is unbounded for ri —)■ 0 while 
the left hand side is bounded above for r* sufficiently small. 

It follows that the diagonal D and the sptP do not intersect and since both sets 
are compact inside Q{r), they must have positive distance in Q{r) □ 


Let Pn denote a sequence of minimisers of K{cn, •) converging to a minimiser P 
of K{c, •). We choose 0 <71 such that P{Q{71)) = M > 0 and we consider a{7l) 
according to Proposition 13.81 above. Since Q{7V) is open we have, 


liminfF„(Q(P)) > M 

n^oo 


and then for n big enough Pn{Q{7Z)) > —. Since Da{Ti) H Q(7l) is closed we have 

lim Pn{Da(Ti) n Q{71)) = 0. 

n—>-CxD 


Then for n big enough 

^ < F„(< 3 (K) \ 


(3.7) 


Proposition 3.10. Assume that p does not have atoms and let n satisfies \3. ?)) 
above, then there exists a maximiser u of D{cn,-) and two positive constants 
and kn such that 


u 


< 


2N{N -If 

Tn 


{N-lfk^. 
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We will later show that we can control r^, and kn nnifornily in n. 


Proof. Consider a point (xi,..., xn) G sptPn \ -Da( 7 ^) and in QiJZ). 

We start from a maximiser (mi, ... ,U]sf) of fl3.6p and we remark that withont 

loss of generality we may assnme that Ui{xi) = ■= kn for all i and 


that 


N 

Ui{x) = inf{c(i/i,..., X, Vi+i ..., ?/ 7 v) - ^ Uk{yk)}-, 


for all i and x. Then we begin by choosing < 


a(7^) 


and we have the following 


estimate. If x ^ U^2-®(^u^n,) then 

Ul{x) < Cn{x,X2, ■ ■ ■ ,Xn) - U2{X2) - Un{xn) < 

< c{x,X2, . . . ,Xn) - U2{X2) - Un{xn) < 

N{N-1) 


< 


{N - l)kn. 


We now select a, possibly smaller, so that 

p{B{xi, r„)) H-h p{B{xn, r„)) < £ << 


M 


(3.8) 


(3.9) 


It follows that the set 

{y G Q{n) n sptPn I y ^ Daiji), Vi ^ B{xj, r„) for i,j = 1,... ,N} 

has positive measnre. Next, take (^i,... in this last set. Since Pi does 
not belong to the balls centered at Xj the estimate above holds and then 

^2(^2) H- = Cn{yi,...,yN)-ui{yi)> 

N{N -1) 


> 


{N-l)kn). 


Finally, np to a division by 2 of we have that for all x G u(l 2 -S(Ti, r„) 

Ul{x) < C„(x,y 2 , ■ ■ ■ ,^jv) - “ 2 (^ 2 )-Itiv(^Ar) < 

< c(x, I/ 2 , • • • , Hw) - “ 2 ( 1 / 2 )-“Y(y 7 v) < 

A^(A^-l) N{N-l) 


< 


-{N-l)kr, 


(3.10) 


This completes the estimate from above of Ui. The same compntation holds for 
the other -Uj. The estimate from above of the Ui given by 03.81) and O3.10p translates 
in an estimates from below which holds p-a.e.. Indeed for p—a.e. x there holds 


ni(x) = inf{Cn(x,X2, . . . ,X„)-M2 (t2)- un{xn)} > {N-l)^kn 


2N{N -ly 
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It remains to remark that 


N 




Z =1 


is a Kantorovich potential for c„ and it satisfies the required estimate. □ 

Proposition 3.11. The constants and kn in Proposition\3A^ can he controlled 
uniformly in n. 

Proof. Consider a point x = (xi,... ,xn) G sptP \ Da{n) and in Q{TZ). Then 
consider r such that 
a(7^) 


(1) r < 


M 


( 2 ) NY,p(B{x.,r) <-. 


Now we will show that in the construction of Proposition 13.101 we may choose t 

c(xi, . . . , X^r) r * 

- — - and = -. Infact, since Pn ^ P we may choose in the estimates 

of Proposition 13.101 a sequence of points x” = (x^,... ,x^) —t (xi,..., xat). The 
convergence of x” to x already gives the required convergence of kn- Moreover, 

T 

we also have that for n big enough B{xf, -) C B{xi,r) for all i and then fl3.9p 

is satisfied. Finally in Proposition 13. lOl we divided by 2 to have some distance 
_ _ r 

between Xj and yj and this last division brings us to = -. □ 

3.2. Conclusion and remarks. 

Theorem 3.12. The following duality holds 

1 


mm / 

n(p) jRJVd 


E 


Xi 


Xi 


-(iP(xi, ...,xn) 


1<2<J<A^ 

= sup I N u{x)dp{x) : u & and u{xi) + • • • + u{xn) < c(xi,..., x„) 


(3.11) 


and the right-hand side of equation above admits a bounded maximizer. 

Proof. By Prop. 13.31 


mm 


E 

1 < 2 <_ 7 <W 


n(p) jRJVd \Xi Xj 


-dP{xi,..., x^) = lim min K'(c„, P). 


By Propositions 13.61 , 13.7113.101 and 13.111 

minK'(c„, P) = maxP(c„, u), 


sup In u{x)dp{x) 


u & Lp and 


r \ . . / ^\ ^ / \ r = lim maxP(cn,M 

u[xi) ^ -h 'U(x„) < c(xi, . . . , Xn) I n^oo 
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and there exists a sequence {un\ where Un is a maximizer of D{cn, u) and where 
|tt„| is uniformly bounded, thus weakly compact in □ 

Corollary 3.13. Assume that p does not have atoms then there exists 0 < a sueh 
that for every minimizer P of K{c, ■) 

P{D^) = 0 

Proof. For every Kantorovich potential u exploiting the duality (or complementary 
slackness) relations we obtains that for every minimizer P of K{c, •) 

u{xi) +-h u{x]y) = c(xi ,... ,xn) P — a.e. 

Since there exists a bounded Kantorovich potential we obtain the conclusion. □ 


Remark 3.14. Driven by the interest in some application we have chosen to present 
the result in the case in the case of Coulomb cost. However hxing some constants 
and exponent (in particular in Lemma 3.4, Prop. 3.7 and Subsection 13.ip and the 
dehnition of bad set the same result may be proved for several costs which are 
lower semi-continuous and bounded from below. Among them 


E 



The case of even less regular costs has been considered in jl] but in that case it is 
necessary to give a different interpretation of the problems. 

In [1^ Kellerer considered a duality theory for multimarginal problems for a 
very wide class of costs. However in order to have existence of maximizers for the 
dual problems (see Theorem 2.21 in El) it is required that the cost c be controlled 
by a direct sum of functions in also from above. In the case of the Coulomb 
cost this would read 

E l- ^ -r < u{xi) H-h u{xisf) 

_ Ay 7 Ay 7 


and this is not possible since the right-hand side allows for xi = ■ ■ ■ = Xn- 
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